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Abstract 

^ , Several distributions are studied, simultaneously in the real, complex, quaternion 

[~ — ' and octonion cases. Specifically, these are the central, nonsingular matricvariate 

OO . and matrix multivariate T and beta type II distributions and the joint density of 

' the singular values are obtained for real normed division algebras. 

in ; ^ 

1 Introduction 

The complex case has renewed interest in multivariate analysis in diverse areas of science and 
technology, see Mehta [l^, Ratnarajah et al. [s^ and Micheas et al. [s^, among many others. 
Moreover, diverse works involving multivariate analysis have appeared in the context of the 
^ ■ quaternion case, see Bhavsar Q, Forrester 17|, Li and Xue 28 1, among others. However, 



H I with respect to the octonion case, only a few, theoretical results have been published, see 

Forrester Il7| and Khatri 25[. This lack of widespread interest may be, because as stated 



by Baez [l|, ...there is still no proof that the octonions are useful for understanding the real 
world. Nevertheless, for the sake of completeness, we include results in the octonion case as 
conjectures. 

Using definitions, properties and notation from abstract algebra, we propose a unified 
approach that enables the simultaneous study of the distribution of a random matrix in the 
real, complex, quaternion and octonion cases, which is generically termed the, distribution 
of a random matrix for real normed division algebras. 

In particular, the matricvariate T distribution has been studied by many authors in the 



real case, see Dickey [12!], Box and Tiao [3], Press |32|, Kotz and Nadarajah [27[ and Diaz- 
Garcia and Gutierrez- Jaimez Q, among many others. The matricvariate T distribution 
appears in the frequentist approach to normal regression as the distribution of the Stu- 
dentised error, see Di'az-Garci'a and Gutierrez- Jaimez (if] and Kotz and Nadarajah 27]. In 
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Bayesian conjugate-prior and diffuse-prior analysis for the same sampling models, it appears 
as the marginal prior or posterior distribution of the unknown coefficients matrix, and also 
as the predictive distribution of a future data array, see Dickey , Box and Tiao , Press 



[32j . Dfaz-Garcfa and Ramos-Quiroga Fang and Li [15| and Kotz and Nadarajah [27 1. 
It has been applied in microeconomic modeling to describe the operation of a market for 
a particular economic commodity and in macroeconomic modeling to describe the interre- 
lations between a large number of macroeconomic variables, as an application of the linear 
simultaneous equation model, see Kotz and Nadarajah (27j . In the complex case, the matric- 
variate T distribution has been apphed in Bayesian estimation of a multivariate regression 
model, see Kotz and Nadarajah [27|. No less important is the role of the T distribution, 
because if the matrix T has a T distribution, then the matrix TT* (or T*T) is distributed 
as beta type II; and the distribution of the latter, in particular, plays a fundamental role in 
the MANOVA model, see Srivastava & Khatri (SJI, Press fs?], and Muirhead [3l|. 

In this work, the nonsingular central matricvariate T and the beta type II distributions 
and some of their basic properties are studied, see Section [51 The matrix multivariate T 
distribution and its corresponding beta type II distribution is studied in Section [3] Finally, 
the joint densities of the singular values are derived in Section [Jj We emphasize that all 
these results are derived for real normed division algebras. Some concepts and the notation 
of abstract algebra and Jacobians are summarised as an Appendix. 



2 Matricvariate T distribution 

We begin this section by distinguishing between matricvariate and matrix multivariate (or 
matrix variate) distributions. We say that the random matrix X has a matricvariate dis- 
tribution if the kernel g{-) of its density is written as a function solely in terms of the 
determinant operator (7(|X|). In any other case it is said that X has a matrix multivariate 
distribution. The term matricvariate distribution was introducing by Dickey [l3|, but the 
expression matrix-variate distribution (or matrix variate distributioii) was later used to de- 
scribe any distribution of a random matrix, see Gupta, and Varga [2l|, Gupta and Nagar 
[20| . and references therein. Alternatively, the term matrix multivariate (instead of matrix 
variate) has been used by Goodall and Mardia [31 , and this is the approach adopted in our 
paper. 

Theorem 2.1. Let T £ C^-^ „ defined as 

where L is any square root ofY such that LL* = V ^ yV/^(i/, H), H g and v > /3{m—l); 
independent ofY ^ A/j^xn(^i ® ^ ^ ^n- Then the density o/ T is 

+ + (T - m)I]-1(T - /x)*r^("+^)/^ (1) 

^™n/3/2r^[^j,/2]|H|/3./2|S|/3™/2' 

which is termed the matricvariate T distributior0 and is denoted as 



Proof. From Kabe [2j| and Diaz-Garcia and Gutierrez- Jaimez [g, [lOj, the joint density of V 
and Y is 

cx |v|/5(''-™+i)/2-i etr{-/3(H-iV + Y,-^Y*Y)/2}, 



^In the literature, it is customary to use the expressions real matricvariate T distribution, complex ma- 
tricvariate T distribution, quaternion matricvariate T distribution and octonion matricvariate T distribution 
distribution; here, however, we use simply matricvariate T distribution as the generic term. 
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where the constant of proportionahty is 

1 



(2/3-i)'9"W2r^[^j,/2]|H|'3W2 (27r/3-i)/5™"/2|S|/3W2 ' 
Making the change of variable Y = LT, where V = LL*, then by (IA-3P 
(dV)(dY) = |LL*|'^"/2(dV)(dT) = |V|'^"/2(dV)(dT). 
Thus, the joint density of V and T is 

cx |v|^'('^+"-"+i)/2-i etr{-/3(H-i + (T - /y)S-i(T - u)*)/2}. 
Finally, integrating over V S we have 

(2^-i)/3m(«+,.)/2pj^(^ + i/)/2]|H-i + (T - iy)S-i(T - ,^)*|-/5("+'^)/2^ 

from which the desired result is obtained. □ □ 

Now, observe that by Dickey 

Ttmn + v)/2] _ Ti[f3{n + v)/2] 

7r™"'5/2r^[/3z//2] 7r™"^/2r^[/3(n + i/ - to)/2] 

and 

|H-i + (T - iy)S-i(T - = |H|-i|S|-i|S + (T - i.)*H(T - 
from which, alternatively, the density ([l} can be expressed as 

rg[/^(n + .)/2]|S|^V2|E|/^(»^---)/2 ^ _ _ 
7r™"'3/2rg[/3(„ + j,_^)/2] 

Corollary 2.1. Lef T e -C^ „ (ie/?nerf as 

T = XLj^i + /X 

w/iere Li is any square root of U such that LiL^ = U ~ V^n{v + n — m, S € 

independent ofX. ^ M^y.j^{0,a~^ In), witha£^f^^. Then, T ^ 7^x„(;/, /x, H, S). 

Proof. The proof is a verbatim copy of the proof of Theorem 12.11 □ □ 

Now, assume that T ^ 7^5^„(z^, 0, Im, In) with n > rn and let F e defined as 
F = XT* then, under the conditions of Theorem 12.11 and Corollary [231 we have 

F = L-^YY*(L-i)* = L-iS(L-i)* 

= xu-^x*, 

with S = YY* - Wi{n,Irn), n > P{m - 1). Thus: 
Theorem 2.2. The density ofF is 

1 ■p|/3(«-m+l)/2-l|j^^^ ^ p|-/3(n+i/)/2^ ^3^ 



Bi[Pi^/2,l3n/2] 



where B^[-, ■] is given by HA-S^) and F is said to have a matricvariate beta type II distribu- 
tion. 
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Proof. The proof follows from ((T)) by applying (|A-ip and (|A-6|) . □ □ 

In addition, assume that n < m and let F G defined as F = T*T then, under the 
conditions of Theorem 12.11 and Corollary 12. II we have 

F = X*U-iX 

^ L-iX*X(L-i)* = L-iSi(L-i)* 

with Si = X*X - W/^(m, I„), m> /3{n~ 1), Thus: 

Theorem 2.3. F has the density 

I |F|/5('"-«+l)/2-l|I _^ p|-/3(n+i.)/2-l ^^-v 

B^[/3{iy + n-m)/2,(3m/2y ' " ' 

Furthermore, we say that F /las a matricvariate beta type II distribution. 



Proof. The proof is the same as that given in Theorem l2.2l Alternatively, observe that den- 
sity dU can be obtained from density ^ making the following substitutions, see Muirhead 



M Eq. (7), p. 455] and Srivastava & Khatri [3J, p. 96], 

TO— >-n, n— )-TO, v ^ v + n — m. □ (5) 

□ 

Finally, assume that M e >Cmxm ^-^Y square root of constant matrix A = MM* e ^P^. 
Also, define Z = M*FM, therefore: 



Corollary 2.2. The density of Z is 



2;|/J(n-m+l)/2-l|^ ^ 2|-'3("+'')/2, 



S^[M2,/3n/2] 

Z is said to have a nonstandardised matricvariate beta type II distribution. 

Proof. The proof follows from ([3]) by applying (|A-4p . □ □ 

Densities ([3|) and (H]) have been studied by several authors in the real case, see Khatri 
26l and Srivastava & Khatri [s^, Gupta and Nagar [i^], among many others; and by James 
23j , Muirhead [3l|, Diaz-Garcia and Gutierrez- Jaimez Q and Dfaz-Garcia and Gutierrez- 
Jaimez Q and Diaz-Garcia and Gutierrez- Jaimez Q, in the noncentral, doubly noncentral, 
singular and nonsingular and complex cases, among many other authors. 

3 Matrix multivariate T distribution 

Theorem 3.1. Let Ti = S-^/^Y+(i e C^^^,^ where {S'^l^f = S^T^{i^,p) (that is, S has a 
gamma distribution with parameters v and p), p > 0, independent of \ ^ J\f^^^^^{0, Im^S), 
S S ^P^. Then the density of Ti is 

rf[/3(z.-KTOn)/2]pW2 1 ...rj. ..-0i.+mny2 

which is termed the matrix multivariate T distribution and is denoted as Ti ^ ■^Tmxni'^^ 
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Proof. The joint density of S and Y is 

cx s^''/2-ietr{-^(s/p + trYY*)/2}, 

where the constant of proportionahty is 

1 1 

c = 



Noting that by (|X3)) 

{ds)idY) = s^"'''/\ds){dTi), 
the desired resuh is obtained analogously to the proof of Theorem 12.11 □ □ 

Corollary 3.1. Assume that Ti ^ A4T^^^{i', fi,I„i,In) , and let M and N be any square 
root of the constant matrices A = MM* £ and A = NN* G respectively. Also let 
(M*)-iTiN-i +n = Qie C^^,„ where n e >C^j „ is constant. Then, 

[1 - tr A(Q, - M)A(Q, - ^)T^(^+™")/^ , 

Hence, we write Qi ^ A^T^xnl'^i /^i -'^)- 
Proof. The proof follows observing that, by (jA-3p 

(dTi) = |MM*|'5"/2|NN*|'5™/2(^Q^) ^ |A|'^"/2|^|/3m/2(^Q^)^ □ 

□ 

Now, assuming that Ti ~ A^7"(^xn(^i 0^ Im, In), with n> m and defining Fi = TiT^; G 
then, under the conditions of Theorem 13. II we have that 

Fi = S-^^YY* = S^^W 

where W = YY* - WiUn, Im), n > p{m - 1). 
Theorem 3.2. The density of Fi is 

-■13 



Tj^lPjv + mn)/2] ,^(„_„+i)/2_w-|^ ^-P{^+mn)/2 

rf[M2]r^[/3n/2]' " " 



(7) 



Fi is said to have a matrix multivariate beta type II distribution. 

Proof The proof follows from © by applying (|A-ip and (|A-6|) . □ □ 

Similarly, if n < m and Fi = T^Ti e ^!^. 
Theorem 3.3. Fi has the density 

rf[M2]r^[/3m/2]' 

Thus, Fi is said /laiJe a matrix multivariate distribution type II distribution. 
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Proof. The proof is the same as that given in Theorem 13.21 Alternatively, the density ([5]) 
can be obtained from density ([7]) by making the following substitutions, 

m — >■ n, n — !■ m. □ (9) 

□ 

As in the matricvariate beta type II distributions, assume that A G ^mxm square 
root of constant matrix 11 = AA* G Also, define Z = A* FA, therefore: 



Corollary 3.2. The density of Z is 

rf[/3(i/ + mn)/2]|np"/2 



z|/3("-™+i)/2-iQ ^_ nz)"'^("+''^/^. 



r^[/3z//2]r^[/?m/2] 

Z is said to have a nonstandardised matrix multivariate beta type II distribution. 

Proof. The proof follows from ([7]) by applying (jA-4p . □ □ 

In the real and singular case, the matricvariate and matrix multivariate T distributions 
have been studied by Dfaz-Garcfa and Gutierrez- Jaimez Q. 

4 Singular value densities 

In this section, the joint densities of the singular values of matrices T, T, Ti and Ti are 
derived. In addition, and as a direct consequence, the joint densities of the eigenvalues of 
F, F, Fi and Fi are obtained for real normed division algebras. 

Theorem 4.1. Let Si,. ..,6m. be the singular values of T ^ T^xnC*^' I")' Si > ■ ■ ■ > 
5m > 0- Then its joint density is 



r^[/3m/2]e^[M2,/3n/2]^iV 



where r is defined in Lemma\7S[ 

Proof. This follows immediately from ([1]) , first using (jA-51) and then applying (jA-ll) . □ □ 



The joint density of the singular values of T is obtained from (|T0|) after making the 
substitutions ([5]). 

Theorem 4.2. Assume that Ti ~ AiT'^^jj^{iy,0,Im,ln) o,nd let ai, ... ,am> ai > ■ ■ ■ > 
, am > 0, be its singular values. Then its joint density is 



2 2n/3 



oti - a 



(11) 



Proof. The proof is identical to that given for Theorem 14.11 □ □ 

Analogously, the joint density of the singular values of Ti is obtained from ([TT|). making 
the substitutions (|9]). 
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Finally, observe that Si — \/eigj(TT*) and ai — ^eigj(TiT*), where eigj(A), i = 
l,...,m, denotes the i-th eigenvalue of A. Let Ai = eigj(TT*) — eigj(F) and ji = 
eigi(TiTJ) = eigj(Fi), observing that, for example, 6i = \f\i. Then 

m mm 

i—1 i—1 i—1 

the corresponding joint density of Ai, . . . , Am, Ai > • • • > \„i > is obtained from (|10p as 

7r -r-r /3(„-m+l)/2-l., , x . w/3(i'+n)/2 TT. v _ ^. 



rf„[/3m/2]Sf„[M2,/3n/2] f^-^ 



n Af + A,)-''(''+")/2 [](A, - A,)^ 



Analogously, the joint density of 71, ... , 7^, 71 > • • • > 7m > 0, is obtained from (jlip as 

-,(3(i' + mn)/2 



Remark 4.1. Observe that Y e n ^ matrix multivariate elliptically contoured 
distribution for real normed division algebras if its density, with respect to the Lebesgue 
measure, is given by (see Diaz-Garcia and Gutierrez- Jaimez [^): 

^'("^'"^ ,/.{tr[E-i(Y-/x)0-i(Y-/xr]}, (12) 



|5]|/3«/2|0|/3m/2 



where /x S £m^„, S e g The function h : ^ [0,oo) is termed the generator 
function, and it is such that J^fi ul^^^~^h{v?)du < 00 and 



27r/3mr!/2 W^^ 

Such a distribution is denoted by Y ~ v,n (Mi 0, /i), for the real case, see Fang and 
Zhang fi^ and Gupta, and Varga [21]; and Micheas et al. [30i| for the complex case. Observe 
that this class of matrix multivariate distributions includes normal, contaminated normal, 
Pearson type II and VII, Kotz, Jensen-Logistic, power exponential and Bessel distributions, 
among others; these distributions have tails that are more or less weighted, and/or present 
a greater or smaller degree of kurtosis than the normal distribution. 

Assume that Y = ( Yi Y2 ) - £,^,^„+^(0, h), n,v > m; and define, 

rnxn mxi' 

T = L^^Yi, where L is any square root of V = Y2Y2 such that LL* — V. Then 
T - 7;fx„(i^,0,I™,I„). From ^ the density of Y is 

(m, n + u)h {tr (Yi Y* + Y2 YJ)} . 

Let V = Y2Y* then by (O, (dYz) = 2-"|V|'5(''-"+i)/2-i((iV)(HidH*). Thus, the 
marginal density of Yi and V is obtained by integrating over Hi e V/^j „ by using (jA-ip . 
hence 

r^[/3i./2] 

Now, let T = L^^Yi, where LL* = V, then by ^ 

(dYi)(dV) = \Vf'l^{dT){d\). 
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Therefore, the joint density of T and V is 



C/^(m, n + z/)7r^-W2 ^^^^j,„ ^ TT*)V} . 

The desired resuh follows by applying the next equally, 

see Di'az-Garcfa and Gutierrez- Jaimez and Diaz-Garcia and Gutierrez- Jaimez [l^]- □ 
Observe that in this case, Yi and Y2 (or V = Y2Y2) are stochastically dependent. 
Furthermore, note that only when the particular matrix multivariate elliptical distribution 
is the matrix multivariate normal distribution, are Yi and Y2 (or V = Y2Y2) indepen- 
dent. Therefore, we can say that the matricvariate T distribution is invariant under the 
family of matrix multivariate elliptical distributions for real normed division algebras, and 
furthermore, its density is the same as when normality is assumed. Analogously, it can be 
proved that the matrix multivariate T, matricvariate and matrix multivariate beta type II 
distributions are invariant under the family of matrix multivariate elliptical distributions 
for real normed division algebras. Furthermore, this invariance prevails under other classes 
of elliptical models for real normed division algebras, see Fang and Zhang [l6j, Gupta, and 



Varga [2l| and Diaz-Garcia and Gutierrez- Jaimez 



Conclusions 

Any topic in statistical literature, is usually first studied in the real case, later in the complex 
case, later for the quaternion case and exceptionally for the octonion case. From the results 
presented in this paper, the real, complex, quaternion and octonion cases are obtaine d by 



simply replacing /3 with 1,2,4 or 8, respectively. Furthermore, as observed by Kabe [24| . 
these results can be extended to hypercomplex cases, that is, for biquaternion and bioctonion 
algebras. Then, from the results presented here, the hypercomplex cases are obtained by 
replacing /3 with 2/3. 
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Appendix 

A detailed discussion of real normed division algebras may be found in Baez f]\ and Gross 
and Richards [l9| . For convenience, we shall introduce some notation, although in general 
we adhere to standard notation forms. 

For our purposes, a vector space is always a finite-dimensional module over the field 
of real numbers. An algebra 5" is a vector space that is equipped with a bilinear map 
m : ^ X ^ ^ ^ termed multiplication and a nonzero element 1 G 'S termed the unit such 
that m(l, a) = m{a, 1) = 1. As usual, we abbreviate m(a, 6) = ab as ab. We do not assume 
^ associative. Given an algebra, we freely think of real numbers as elements of this algebra 
via the map w i— > luI. 

An algebra 5 is a division algebra if given a,b G d with ab = 0, then either a = or 
b = 0. Equivalently, 5^ is a division algebra if the operation of left and right multiplications 
by any nonzero element is invertible. A normed division algebra is an algebra ^ that is 
also a normed vector space with \\ab\\ = ||a||||&||. This implies that 5^ is a division algebra 
and that 1111] = 1. 

There are exactly four normed division algebras: real numbers (iJi), complex numbers 
(£), quaternions (io) and octonions (D), see Baez [ij. We take into account that SR, €, and 
D are the only normed division algebras; moreover, they are the only alternative division 
algebras, and all division algebras have a real dimension of 1, 2, 4 or 8, which is denoted by 
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/3, see Baez [l], Theorems 1, 2 and 3]. In other branches of mathematics, the parameters 
a = 2/ j3 and t = /?/4 are used, see Edelman and Rao [31 and Kabe 24|, respectively. 

Let £^ „ be the hnear space of all m x n matrices of rank m < n over ^ with m distinct 
positive singular values, where ^ denotes a real finite- dimensional normed division algebra. 
Let 5™^" be the set of all m x n matrices over ^. The dimension of g^™^" over jR is /Smn. 

T 

Let A g g^™x"^ then A* = A denotes the usual conjugate transpose. 

Table [T] sets out the equivalence between the same concepts in the four normed division 
algebras. 



Table 1: Notation 



Real 



Complex 



Quaternion 



Octonion 



Generic 
notation 



Semi-orthogonal Semi-unitary Semi-symplectic 

Orthogonal Unitary Symplcctic 

„ , . . Quaternion 

Symmetric Hermitian 7 ■ ■ 

nermitian 



Semi-exceptional 
type 
Exceptional 
type 
Octonion 
hermitian 



U'3(m) 



6: 



/3 



In addition, let be the cone of positive definite matrices S e ^"my-m.^ then is an 
open subset of 

Let 2)^ be the diagonal subgroup of ,„ consisting of all D e g^mx™^ = diag(di, . . . , dm). 
For any matrix X g ^"x™^ denotes the matrix of differentials (dxij). Finally, we 
define the measure or volume element (dX) when X G 5^™^", 6^^, 2)^ or V£ „, see Dimitriu 

If X € g-™xn ^jjgjj (dX) (the Lebesgue measure in ^•'"xn-j ^jgnotes the exterior product 
of the /3mn functionally independent variables 



(dX) = /\ /\ dxij where dxij = ^ dx 
1=1 j=i 



(k) 



fe=i 



If S e 6f„ (or S e 1^(m) is a lower triangular matrix) then (dS) (the Lebesgue measure 
in 6^ or in T^(m), the set of lower triangular matrices) denotes the exterior product of 
the m(m + l)/3/2 functionally independent variables (or denotes the exterior product of the 
m(m — l)/3/2 -|- n functionally independent variables, if su G 3? for alH = 1, . . . , m) 



(dS) 



m /3 

A A^^ 

i<j k=l 



(fc) 
ij 1 



1=1 



The context generally establishes the conditions on the elements of S, that is, if Sij e 9?, 
£ £, £ i3 or g D. It is considered that 



m ^ m m 13 

= A w... A A ^sg^ 

i<j k—1 i—1 i<j k—1 



{dS) = /\ /\dsif^ /\dsu/\ /\ds[ 



Observe, too, that for the Lebesgue measure (dS) defined thus, it is required that S g ^ 
that is, S must be a non singular Hermitian matrix (Hermitian definite positive matrix). 
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If A G 2)^ then (dA) (the Legesgue measure in J)^) denotes the exterior product of the 
l3m functionally independent variables 

71 /3 

(dA) = /\/\dAf' . 

If Hi e V^,_„ then 

m n 

(HJdHi) = /\ /\ h*dh,. 

where H = (HjIKj)* (hi, . . . , hm|hm+i, . . . , h„)* G il''(n). It can be proved that this 
differential form does not depend on the choice of the H2 matrix. When n = 1; 1 defines 
the unit sphere in ^J™. This is, of course, an (m — 1)13- dimensional surface in ^J™. When 
n = m and denoting Hi by H, (HdH*) is termed the Haar measure on il^(m). 
The surface area or volume of the Sticfcl manifold V/^ „ is 

Vol(V/^,J ^ / (HidHD = (A-1) 

where r^[a] denotes the multivariate Gamma function for the space 6^, and is defined by 



f etr{-A}|A|"-("-i)'3/2-i(dA) 

m 



1=1 



where etr(-) = exp(tr(-)), | • | denotes the determinant and Re(a) > (m — l)/3/2, see Gross 
and Richards [l^l ■ Similarly, from Herz the multivariate beta function for the space 
can be defined as 

Bi[b,a] = / |Br-(™-i)''/2-i|I^-B|''-(™+iW2-i((iB) 

Jo<B<I„, 

f |A|'^-(™-i)''/2-i|i^ A|-('^+^)(dA) 

(A-2) 



r5^,.[a + 6] 

where A = (I - B)-^ - I, Re(a) > (m - l)/3/2 and Re(&) > (m - l)/3/2. 

Now, we show two Jacobians in terms of the /3 parameter, which are based on the work 
of Kabe [ij] and Dimitriu ^.S^These results are proposed as extensions of real, complex or 



quaternion cases, see James [23|, Khatri 25|, Mehta [29[, Ratnarajah et al. 33| and Li and 
Xue [2^. 

Lemma .1. Let X and Y G £(^_„; and let Y = AXB + C, where A G /^f„^m, B G 
and C G £fL „ are constant matrices. Then 



(rfY) = |A*A|'^"/2|B*B|^™/2(^x). (A-3) 

V 

constant matrices. Then 



Lemma .2. Lei X and Y G ^3^, and let Y AXA* + C, where A and C e „j are 



(dY) = |A*A|'3("-i)/2+i(dX). (A-4) 
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Lemma .3 (Singular value decomposition, SVD). Let X G „ , such that X = ViDW* 
with Vi e V^_„, W e il''(rn) and D = diag(di, • • • , dm) G 2)^^ di > • • • > c^m > 0. T/ien 

(dX) = 2-™7r^ JJdf("-'"+^^-^ Y[{d^i - djf{dB){VidVl){WdW*), (A-5) 

w/iere 

0, /3 = 1; 
^1 -m, 13 = 2; 
] -2rn, „fl = 4; 
-4m, ^ = 8. 

Lemma .4. Let X e £^_„, and S = XX* e Then 

(rfX) = 2-'"|S|^("-™+i'/2-i(dS)(VidVn. (A-6) 
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